The broad Fe Ka line and supermassive black holes 

Predrag Jovanovic 3 '* 

"Astronomical Observatory, Volgina 7, 1 1060 Belgrade, Serbia 



Abstract 



o 

O 



Here we present an overview of some of the most significant observational and theoretical studies of the broad Fe Ka spectral line, 
which is believed to originate from the innermost regions of relativistic accretion disks around central supermassive black holes 
of galaxies. The most important results of our investigations in this field are also listed. All these investigations indicate that the 
broad Fe Ka line is a powerful tool for studying the properties of the supermassive black holes (such as their masses and spins), 
space-time geometry (metric) in their vicinity, their accretion physics, probing the effects of their strong gravitational fields, and for 



testing the certain predictions of General Relativity. 
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1. Supermassive black holes 
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According t o General Rela tivity (Einstein, 19161) . a black 
hole (named by Wheeler , 19681) is a region of space-time around 
a collapsed mass with a gravitational field that has became 
so strong that nothing (including electromagnetic radiation) 
can escape from its attraction, after crossing its event hori- 
zon (WaldJ, 119841) . The problem of s uch catastrophic gravi- 
tational collapse was first adressed by IChandrasekhari (1193 II) 
when he discovered the upper mass limit for ideal white dwarfs, 
comp osed of a degenerate electron-gas. Subseq uen t studies 
by|Landauldl932h : IOppenheimer & Snvderid 19391) and lPenrose 
(1965) established the modern theory of gravitational collapse. 
The theory of black holes, as well as their fundamental prop- 
erties are presented in numerous mono graphs and papers (see 
e.g. Cartel, 1973 : Chandrasekharl 1983b . and therefore in this 
review we will focus only on those studies which are necessary 
for the discussion of the relation between black holes and the 
Fe Ka spectral line. 

All black holes in nature are commonly classified accord- 
ing to their masses as: supermassive black holes (with masses 
M BH ~ 10 5 - 10 10 M ), intermediate-mass black holes (M BH ~ 
10 2 - 10 s M ), stellar-mass black holes (M BH < 10 2 M G ), 
mini and micro black holes (Mbh ^ M B ). A crucial event 
for the acceptance of bl ack holes was the discovery of pulsars 
by iHewish et al.1 dl968l) . because it was the clear evidence of 
the existence of neutron stars, and therefore, confirmation of 
Chandrasekhar limit. The first detection of a solar mass black 
hole came in 1972, when the mass of the rapidly variable X-ray 
source Cygnus X-l was proven t o be above the maximum al- 
lowed for a neutron star (see e.g. iFerrarese & Fordl 2005, and 
references therein). 



Nowadays, it is widely accepted that supermassive black 
holes are located in the centers of most of galaxies, and thus 
have a fundamental influence on galactic formation and evo- 
lution. Some of the first indirect arguments that they ex- 



lst in galactic nucle i are g 
Lvnden-Bell & Reesl d 1 97 lb . 



ley 

ist in galactic nucle i are given by iLvnden-Belll ([1969) and 



Also, accordin g to the unifica 



tion model of active galactic n uclei - AGN dAntonuccil 11993 



Petersonl 1 19971: iKroliki [l999), they are most likely powered 
by the accretion of gas onto their central s upermassive black 
holes with mass rang ing from 10 5 to 10 9 M G ( Kaspi et al. . 2000; 



Peterson et al., 2004). 



1.1. Space-time geometry in vicinity of supermassive black 
holes 

In general, black holes have three measurable parameters 
(not including the Hawking temperature): charge, mass (and 
hence gravitational field) and angular momentum (or spin). In 
the case of supermassive black holes, only the latter two are of 
sufficient importance because they are responsible for several 
effects which can be observationally detected (for more details 
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see e.g. Uovanovic & Pop ovic, 2008 

A central supermassive black hole of a galaxy can be sta- 
tionary or rotating. In the first case, space-time geometry in the 
black hole vicinity depends on ly on its mass M and is described 
by the Schwarzschild metric ( Schwarzschild , 19161) . while in 
the latter case, space-time geometry depends also on spin a (i.e. 
angular momentum J = aM) of the black hole and is described 
by Kerr metric (IKernll963l) . Therefore, Schwarzschild metric 
describes the spherically symmetric gravitational field around a 
time-steady non-rotating black hole in vacuum, while the Kerr 
metric describes the gravitational field outside an uncharged 
rotating black hole and, contrary to the Schwarzschild metric , 
is no longer spherically symmetric. iBover & Lindquist] d 19671) 
discovered a generalized coordinate syste m in which Ke rr met- 
ric is most commonly used today (see e.g. lCarten,ll973l) . 
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Several characteristic radii can be defined around 
black holes, and t he m ost important are (see e.g. 
Jovanovic & Popovic , l2009t) : 1 . Schwarzschild radius 
Rs = 2GM/c 2 (representing the limiting radius below which 
a collapsed mass form a spherically symmetric non-rotating 
black hole), where G is the gravitational constant and c is the 
speed of light; 2. gravitational radius R s (being a half of Rs 
and usually used as a unit for distance around a black hole); 3. 
radius of event horizon Rf, (representing space-time boundary 
below which events cannot affect an outside observer); and 
4. radius of marginally stable orbit R ms (representing the 
minimum allowed radius of a stable circular equatorial orbit 
around a black hole). 

The equations g overning phot o n trajectories in the K err met- 



ric are derived by ICarter (1968). Misner et al. ( 1973b showed 



that certain relativistic effects in vicinity of a black hole, such 
as light banding, gravitational and Doppler shifts as seen by dis- 
tant observer could be computed by solving the equ ation of the 
photon orbit. ICunningham & Bardeenl (119721 1 1973b calculated 
radiation of a point source in a circular orbit in the equatorial 
plane around an extreme Kerr metric black hole (see Fig. [TJ. 




Figure 1 : A circular orbit (solid line) of a point source in the equatorial plane 
around an extreme Kerr metric black hole. Image Credit: Cunningham & 
Bardeen, Copyright: ApJ, 183, 237 (1973). 



Bardeen (1970) studied the effects of accretion from a disk 



of gas orbiting a black hole on its properties, and showed that 
the Kerr metric is more appropriate for describing supermassive 
black holes than the Schwarzschild metric. This conclusion is 



strengthened even more by ThorngJ (11974b . who showed that a 
black hole at the center of an accretion disk would be spun up to 
the maximum possible value of a « 0.998 in a relatively short 
time. In such a case the disk would be extended down to about 
1.23 R g , while in the case of a non-rotating black hole it could 
extend only down to 6 R g . 

1.2. Accretion onto the supermassive black holes 

Some obs erved quasars have luminosities of up to Lb \ ~ 
10 47 erg/s ( Sch neider , 120061) . The corresponding total energy 
emitted during the lifetime of such quasars can be estimated 
to E > 3 x 10 61 erg, assuming that their minimum age is about 
10 7 years and tha t their luminositie s do not change significantly 
over the lifetime dSchneidenl2006l) . Accretion is the only mech- 
anism which can yield sufficient efficiency e (defined as the 
mass fraction of fuel that is converted into energy, according 



to E = emc 2 ) to explain such high energies. The maximum ef- 
ficiency of accretion is about e ~ 6% for a non-rota ting black 
hole and e ~ 29% for a maximally rotating one ([Schneider, 
2006b . 



On the other hand, if the cold matter which was initially at 
rest and without magnetic field, was subjected to free radial in- 
fall it would accrete to the central black hole without any energ y 



release or observational effects (Sh akura & Sunvaevi 11973b . 
However, in the case of AGN the accreting matter has a signifi- 
cant angular momentum which does not allow its free infall. At 
the marginally stable orbit of central black hole with mass of 
10 8 M o , the specific angular momentum is * 1 x 10 24 cm 2 /s, 
which is much less even in comparison to a typical galaxy 
where the speciffic angular mom entum of orbiting material is 
« 6 x 10 28 cm 2 /s dKrolikl Il999h . It means that approach of 
accreting material toward a black hole requires a loss of the 
greatest fraction of its initial angular momentum. Mechanisms 
which can contribute to such loss of angular momentum are vis- 
cosity, nonaxisymmetric gravitati onal forces, magnetic forces 
( Krolik , 1999)) . outflowing winds (Beg elman et al. , 1983 : Elvis , 
2000b . or the Magnetorotation al I nstabihty (MRI) which was 



discovered by Velikhov ( 1959) and Chandrasekhar ( 1960), but 



did not come to the attention of the as t rophysical c ommunity 
until rediscovered bv lBalbus & Hawlevl (1 1 99 lL 1 1 998h . 



Since the orbit of minimum energy for fixed angular momen- 
tum in any spherically symmetric potential is a circle, the infall 
of accreting material due to loss of its angular momentum will 
be in form of successively smaller and smaller concentric cir- 
cles. Matter traveling along orbits inclined to each other will 
eventually collide in the plane of intersection, and as a result, 
the angular momenta of different gas steams will be mixed and 
equalized. Consequently, all accreting matter will orbit in a 
single plane and will have the same specific angular momen- 
tum at any given radius, meaning that accretion is m ost likely 
performed through an accretion disk ([Krolik, Il999t) . The as- 
sumption of a disk geometry for the distribution of the X-ray 
emitters in the central parts of AGN is al so supported by the 



spectral sha pe of the observed Fe Ka line (Nand ra et al 
Il999ll2007b. 



1997 



Shak ura & Sunvaevi d 1973b developed what is now called 
"standard model of accretion disk", in which accretion occurs 
via an optically thick and geometrically thin disk and where 
the spectrum of thermal radiation emitted from the disk surface 
depends on its structure and temperature, and therefore on the 
distance to the central black hole. The distribution of the tem- 
perature along the radius of the disk i n this model is presented 
in Fig. [2] (see Popovic et al. , 2006al for more details). This 
model was originally developed to describe the accretion disks 
around stellar sized black holes in the binary systems, but with 
certain modifications it could be also applied on accretion disks 
around central supermassive black holes of AGN. 

The total energy release of emitted radiation is mainly de- 
termined by the rate of matter inflow into the disk on its outer 
boundary, i.e. by its accretion rate M. If the accretion converts 
matter to radiation with fixed radiative efficiency rj (in rest-mass 
units), then a characteristic scale for the accretion rate is the Ed- 
dington accretion rate for which the total release of energy in 
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Figure 2: The distribution of the temperature as a function of the disk radius R, given for two different values of angular momentum a. Negative values of R 
correspond to the approaching and positive values to the receding side of the disk. Image Credit: Popovic et al., Copyright: ApJ, 637, 620 (2006). 
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Figure 3: Magnetohydrodynamic simulations of an accretion disk as seen by a distant observer at the following inclination angles I Armitage and Reynolds, 2003): 
5° (upper left), 30° (upper right), 55° (lower left) and 80° (lower right). The disk rotation is in counterclockwise direction. Image Credit: Armitage & Reynolds, 
Copyright: MNRAS, 341, 1041 (2003). 



the disk L = rjMc 2 is equal to the Eddington luminosity, a crit- 
ical luminosity beyond whi ch the radiat i on force overpowers 
gravi ty for any given mass dKrolikl 1 1999c IShakura & Sunvaevt 
1973b . Observed AGN have luminosities from ~ 10 42 to ~ 10 48 
erg/s, which means that the ir central bla ck holes must have 
masses from 10 5 to 10 9 M Q dKrolikl 1 1999h . 

Using the Faint Object Spectrograph (FOS) on the Hubble 



Space Telescope (HST), Har ms et al.1 (1994) found strong evi- 
dence for a disk of ionized gas around a supermassive nuclear 
black hole in M87. 



lArmitage and Reynolds ( 2003 ) performed magnetohydrody- 
namic simulations of an accretion disk as seen by a distant ob- 
server at four different inclinations (see Fig. [3). These simula- 
tions showed the existence of especially bright regions in form 
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of arcs within the turbulent flow, which trace out the photon 
trajectories close to the radius of marginally stable orbit. These 
arcs are much brighter on approaching side of the disk due to 
higher surface temperature, as shown in Fig. [3] and also due to 
Doppler boosting and relativistic beaming which both enhance 
the flux observed in the direction of motion and diminish the 
flux in the opposite direction. The influence of these relativistic 
effects on the shape of the Fe Ka line emitted from the accretion 
disk will be explained in more details in §2.3. 

2. Fe Ka spectral line 

2.1. Observational studies 

Some observational studies suggest that the broad fluores- 
cent/recombination iron Ka line at 6.4 keV originates from the 
innermost part of accretion disk, close to the central black hole. 
For example, using XMM- Newton observation of the broad- 
line radio galaxy 4C 74.26. iBallantvne & Fabianl d2005l) found 
a broad Fe Ka line emitted from a region which had inner ra- 
dius close to the innermost stable circular orbit (R ms ) for a maxi- 
mally spinning black hole and outer radius within 10 R g . There- 
fore, the broad Fe Ka line is an important indicator of accreting 
flows around supermassive black holes. At the same time, it is 
the strongest line of the X-ray radiation, and it can be found in 
the spectra of all types of accreting sources: binary black hole 
and neutron star systems, cata clysmic variable stars a nd AGN. 

It was first discovered by iMushotzkv et a i1 dl978h in OSO- 
8 X-ray observations of Cen A (NGC 5128) during 1975 and 
1976 (see Fig. |4|. NGC 5128 is one of the closest radio galax- 
ies and one of the first extragalactic objects to be identified as 
an X-ray source. The reflection of the X-rays in K lines of 
heavy elements from a cold surface was a l so di scussed in the 
case of the X-ray binaries (see e.g. iBaskoi 1978 ). and Exosat 
Observatory detected a broad emission iron K line at 6.2 keV in 
the X-ray spectrum of Cyg X-l, a binary X-ray sourc e whic h 
was the best known black hole candidate dBarr et all |j985). 
Ginga satellite detected a strong iron K a fluorescence line in 
a number of Seyfert 1 galaxies (see e.g . Kuniedaet all 1990; 



CEN-A EDGE AND LINE 



Matsuoka et al., 1990; Pi ro et all Il990t iPounds et all 11990 



Nandra et al.LI 1 99 ll )7The Fe Ka line of active galactic nuclei, as 
well as different geometries and astrophysical conditions nec- 
essary fo r its emerg e nce, was intensively studied aro und 1990s 



(see e.g. Ilvlatt et all 1 199 lb iGeorge & Fabianl Il99ll and refer 



ences therein). 

The CCD detectors on Japanese ASCA satellite were the first 
instruments with sufficie nt spectral resolution and sensitivity in 
the X-ray band, by which Tana ka et al. I d 19951) obtained the first 
convincing proof for the existence of a relativistically broad- 
ened Fe Ka line in AGN spectra. This discovery was made 
after four-day observations of Seyfert 1 galaxy MCG-6-30-15 
(see Fig. [5). 

Several studies h ave been per f ormed over samples of lo 
cal AGN (se e e.g. iNandraet all Il997t " 



Yaaoobetal., 2005; 



Nandr a et all 120071 iBianchi et all [2008; Markowi tz et al 
20081 etc), as well as from distant quasars dCorral et a" 



in order to characterize the Fe Ka emission. Nandra et al 
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Figure 4: The Fe Ko- line of Cen A (NGC 5128) observed by OSO-8 dur- 
ing 1975 and 1976 (Mushotz kv et all 1 19781) . Image Credit: Mushotzky et al., 
Copyright: ApJ, 220, 790 (1978). 
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Figure 5: The profile of the Fe K g line from Seyfert 1 galaxy MCG-6-30- 
15 observed by ASCA SIS detector I Tanaka et all 1 19951) and the best fit (blue 
solid line) obtained by a model of the accretion d isk in Schwarzschild metric, 
extending between 3 and 10 Schwarzschild radii <Fabianetallll989l) . Image 
Credit: Tanaka et al., Copyright: Nature, 375, 659 (1995). Image generated by 
Dr. Paul Nandra NASA/GSFC. 



(2007) performed a spectral analysis of a sample of 26 type 1 
to 1.9 Seyferts galaxies (z < 0.05) observed by XMM-Newton. 
They found that a relativistic line is significantly detected in a 
half of their sample (54+10%) with a mean equivalent width 
(EW) of ~ 80 eV, but around 30% of selected AGN showed a 
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relativistic broad line that can be explained by the emission of 
an accretion disk. 

Recent Chandra observations revealed a narrow compo- 
nent of t he Fe Kg line in the X-ray spectra of many AGN 
(see e.g. lYaaoob et all |200U Yaaoob & Pa dmanabhan, 2004 



Page et all BOOfl IShu et all 12010). Although the origin of this 



narrow component is still poorly understood, it is thought to be 
produced by emission from material much further from the cen- 
tral black hole, probably in the outermost regions of the accre- 
tion disk, in the Broad Line Region (BLR) or torus. Therefore, 
in the further text the focus will be only on the broad component 
of the Fe Ka line, which originates from the innermost parts of 
an relativistic accretion disk. 

If the Fe Ka line originated from an arbitrary radius of a non- 
relativistic (Keplerian) accretion disk it would have a symmetri- 
cal profile (due to Doppler effect) with two peaks: a "blue" one 
which is produced by emitting material from the approaching 
side of the disk in respect to an observer, and a "red" one which 
corresponds to emitting material from the receding side of the 
disk. Broadening of the Fe Ka line arises mostly from Gen- 
eral Relativistic effects and large rotational velocities of mate- 
rial emit ting near the black hole. Using ASCA satellite obser- 



vations, Nand ra et alj (U.997) found that, in case of 14 Seyfert 



1 galaxies, Full-Widths at Half-Maximum (FWHM) of their Fe 
Ka lines correspond to velocities of « 50, 000 km/s. In some 
cases (e.g. for Seyfert 1 galaxy MCG-6-3 0-15), FWHM veloc - 



ity reaches 30% of speed of light (see e.g. lNandra et a l., 2007) 



It means that in the vicinity of the central black hole, orbital ve- 
locities of the emitting material are relativistic, causing the en- 
hancement of the Fe Ka line "blue" peak in regard to its "red" 
peak (relativistic beaming). Taking into account the integral 
emission in the line over all radii of accretion disk, one can 
obtain the line with asymmetrical and highly broadened profile 
(Fab ian et alll2000l) . The "blue" peak is then very narrow and 
bright, while the "red" one is wider and much fainter (see Fig. 
|3). Besides, the gravitational redshift causes further deforma- 
tions of the Fe Ka line profile by smearing the "blue" emission 
into "red" one (see §2.3 for more details). Since the observed 
Fe Ka line profiles are strongly affected by such relativistic ef- 
fects, they represent a fundamental tool for investigating the 
plasma conditions and the space-time geometry in the vicinity 
of the supermassive black holes of AGN. 

2.2. Production of the Fe Ka line 

A substantial amount of X-ray radiation of AGN is thought to 
be emitted from the hot corona sandwiching the accretion disk 
(see Fig. [6]). The Fe Ka line is produced when the hot corona 
irradiates the relatively cold accretion disk by the hard X-ray 
power law continuum, causing among the rest, photoelectric 



absorption followed by fluorescent line emission (Fab ian et al 
2000). When plasma is subjected to the influence of the hard 



X-ray radiation, one of the two /f-shell (n = 1, where n is the 
principal quantum number) electrons of an iron atom (or ion) 
is ejected following the photoelectric absorption of an X-ray 



(Fabi an et al.L 120001). The thresho ld for the absorption by neu- 
tral iron is 7.1 keV dFabian et al. , 2000). The resulting excited 



state decays when an L-shell (n = 2) electron drops into the K- 
shell, releasing 6.4 keV of energy. This energy is either emit- 
ted as an emission-line photon (34% probability) or internally 
absorbed by another electron (66% probability) which is conse- 
quently ejected from the iron ion (Auger effect). 
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Figure 6: Schematic illustration of inner region of accretion disk around black 
hole, where the Fe Ka line is produced. Image Credit: Wilms et al., Copyright: 
MemSAI, 75, 519(2004). 

The fluorescent yield (i.e. the probability that photoelectric 
absorption is followed by fluorescent line emission rather than 
the Auger effect) is a weak function of the ionization state from 
neutral iron (Fe I) up to Fe XXIII dFabian et all I2000T) . For 
lithium-like iron (Fe XXIV) through to hydrogen-like iron (Fe 
XXVI), the lack of at least two electrons in the L-shell means 
that the Auger effect cannot occur. For He and H-like iron ions, 
the line is produced by the capture of free electrons (recombina- 
tion) and the equivalent fl uorescent yield is h igh and it depends 
on the plasma conditions (IFabian et all 120001) . 

For the neutral iron, the Fe Ka line energy is 6 .4 keV (more 



precis ely, there are two components of the line ( Fabi an et al 



2000): Fe Kai at 6.404 and Fe Ka 2 at 6.391 keV), while in the 
case of ionization, the energy of both the photoelectric thresh- 
old and the Fe Ka line are slightly increased. Even for such 
high ionization states of He and H-like iron ions, the Fe Ka 
line energy is increased only to 6.7 and 6.9 keV, respectively 
(Krolik, 1999). Importance of ionization, X-ray scattering and 
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Compton reflection o n the proper t ies of the observed Fe Ka 
line, was discussed bv lTurner et al.ld 1997b . based upon the anal- 
ysis of ASCA observations of 25 Seyfert 2 galaxies. 

Fe Ka line is pretty narrow in itself, but in case when it orig- 
inates from a relativistically rotating accretion disk of AGN it 
becomes wider due to kinematic effects, and also its shape (or 
profile) is changed due to Doppler boosting and gravitational 
redshift (see Fig. |7). Such broadening of the line is very often 
observed in spectra of Seyfert galaxies and is one of the main 
evidences for the existence of a relativistic accretion disk which 
extends deepl y in the gravitational field of the central black hole 
dZvckiLl2004b . 



2.3. Theoretical studies 

In general, there are several approaches to numerically eval- 
uate the li ne profiles emitted by the accretion disk around 
black hole (Reynolds & Nowak, 2003 [): i) analytically in the 
weak field limit dChen et dlll989tlChen & Halperilll989b and 
in the Schwarzschild metric (Fabian et al., 1989; Stella, 1990; 



Matt et all 119921) . ii) by brute-force using the dire ct integration 



of the photon trajectory in the Kerr metric (e.g. iKaras et al 



1992), iii) by expressing observ ed flux over transfer -function 
containing all relativistic effects ( Cunn ingham. 1975|). and iv) 
using so called ray-tracing method in Ke r r metr i c feao et all 
1994 : iBromlev et ail 1 19971: iFanton et all 1 19971: ICadez et all 
1998). 
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Figure 7: The influence of Doppler and transverse-Doppler shifts, relativistic 
beaming and gravitational redshifting on profile of the Fe Kg- line emitted by 
an accretion disk around a black hole. Image Credit: Fabian et al., Copyright: 
PASP, 112, 1145 (2000). 



Chen et al.1 dl989l) and lChen & Halpernl dl989l) developed an 



emitted from a Keplerian accretion disk in which the relativistic 
effects are included in the weak-field limit. Rigorous relativistic 
calculations of the Fe Ka line profile emitted from a geometri- 
cally thin accretion disk around a stationary black hole, i.e. in 
the case of Schwar z schild metric h ave been carried o ut by e.g. 
Fabian et al.1 dl989h : IStellal dl990b : iMatt et al.1 ri 1992b . The rel- 
ativistic effects on the continuum emission from an accretion 
disk around a rota ting black hole (i.e . for Kerr metric) were 
first calculated by Cunningham ( 1975b. while t he line profiles 
in suc h a case were simulated bv lLaorl dl99lb . Fabian et al 



(2000) studied the influence of Doppler effect and gravitational 
redshift on the shape of Fe Ka line emitted by an accretion disk 
around a black hole and found that, as a result, the line profile 
appears to be broadened and skewed (see Fig. |7). 

The emission from accretion disk can be analyzed also by nu- 
meri cal simulations based on ray-tracing method in Kerr met- 
ric dBao et a 111 994 IBromlev et all 1 1 9971; IFanton et all [1997 



Cadez et al. , 1998), taking into account only photon trajectories 



reaching the observer's sky plane. In this method one divides 
the image of the disk on the observer's sky into a number of 
small elements (pixels). For each pixel, the photon trajectory is 
traced backward from the observer by following the geodesies 
in a Kerr space-time, until it crosses the plane of the disk (see 
Fig. [8}. Then, the flux density of the radiation emitted by the 
disk at that point, as well as the redshift factor of the photon are 
calculated. In that way, one can obtain the color images of the 
accretion disk which a distant observer would see by a high res- 
olution telescope. The simulated line profiles can be calculated 
taking into account the intensities and received photon ener- 



riesof all pixels of the corresponding disk image. Cad ez et al 



1998b developed a variant of ray-tracing method based on the 
pseudo-analytical integration of the geodesic equations which 
describe the photon trajectories in the Kerr metric. 

The complex profile of the Fe Ka line depends on different 
parameters of the a ccretion disk and central black hole (see e.g. 
Fabia n et al. , 2000). In order to study the size of the Fe Ka 



approach for calculating the intensities and profiles of the lines 



line emitting region, as well as its location in the disk, one can 
assume that the line is e mitted from a region in form of a nar- 
row ring. For example, [jovanovi c & Porjovic (2008) assumed 
a line emitting region with width equals to 1 R g , located be- 
tween: a) R in = 6 R g and R ou , = 7 R g and b) R in = 50 R g and 
Rout = 51 R g . These two cases are presented in Fig. [9] From 
Fig. [9] one can see how the Fe Ka line profile is changing as 
the function of distance from central black hole. When the line 
emitters are located at the lower radii of the disk, i.e. closer to 
the central black hole, they rotate faster and the line is broader 
and more asymmetric (see Fig. [9] top-right). If the line emis- 
sion is originating at larger distances from the black hole, its 
emitting material is rotating slower and therefore the line be- 
comes narrower and more symmetric (see Fig. |9]bottom-right). 
In majority of AGN, where the broad Fe Ka line is observed, 
its profile is more similar to the modeled profile as obtained 
under assumption that the line emitters are located close to 
the central black hole dTanaka et al. . 1995 : Nandra et al. . 2007 : 
I Jovanovic & Popovid.l2008b . 

Angular momentum or spin of the central supermassive black 
hole of AGN is a property of the space-time metric and dif- 
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Figure 8: Schematic illustration of the ray-tracing method in the Kerr metric, showing a light ray emitted from some radius of accretion disk in coordinate system 
defined by a rotating black hole with angular momentum a, and observed at a pixel with coordinates (impact parameters) a,fl on the disk image in the observer's 
reference frame. Image Credit: Jovanovic & Popovic, Copyright: Black Holes and Galaxy Formation, Nova Science Publishers Inc, Hauppauge NY, USA, 249-294 
(2009). 




X (in Rg) 



9= E obs/ E ,r 



Figure 9: Left: illustrations of the Fe Ka line emitting region in form of narrow ring with width = lR g , extending from: Rj„ = 6 R ; to R oll , = 7 R g (top) and R,„ = 50 
R g to R OM = 51 R e (bottom). Right: the corresponding Fe Ka line profiles. Image Credit: Jovanovic & Popovic, Copyright: Fortschr. Phys., 56, 456 (2008). 
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X(in^ E ) g = E ohs IE ea 

Figure 10: Illustrations of an accretion disk (left panels) and the corresponding Fe Ka line profiles (right panels) for disk inclination r' = 40° in the Kerr metric for 
angular momentum a = 0. 1 (top) and a = 0.998 (bottom). Image Credit: Jovanovic et al., Copyright: Baltic Astronomy (201 1). 



ferent techniques for its estim ation are developed (see e.g. 
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Figure 1 1 : Modeled Fe Ka spectral line profiles for several values of angu- 
lar momentum parameters a, and for inclination angle i = 40°. Image Credit: 
Jovanovic et al., Copyright: Baltic Astronomy (2011). 



Czernvet all 12011: Daly, 2011). To study the effects of black 



hole spin on the shape of the Fe Ka line we performed sim- 
ulations of disk emission based on ray-tracing method in Kerr 
metric. Simulated images of an accretion disk around a super- 
massive black hole and the corresponding Fe Ka line profiles 
for two differen t value s of its spin are presented in Fig. [ToKsee 



Jovanovic et al.L 1201 ll for more details). In these simulations, 
the disk extends from the radius of marginally stable orbit (R ms ) 
to 40 R g , where R g is the gravitational radius. As one can see 
from Fig. \W\ R ms strongly depends on black hole spin, which 
consequently also affects the profile of the fine originating in the 
inner most parts of the disk (see also iBrenneman & Revnoldsl 



2006). Simulated Fe Ka line profile s for several values of b lack 



hole spin are presented in Fig. [TTI (J ovanovic et al. , 2011 ). It 
can be seen from Fig. QT|that the black hole spin especially 
affects the red wing of the Fe Ka line, and that this wing 
is more exte nded towards lower energies for higher values of 
the spin (see iRevnolds & Nqwaki 120031 : Ijovanovic & Popovid 



2008; Jovanovic et al. , 201 1|). At the same time, the line be- 



comes wider and its red peak brighter. 
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Figure 12: Time-averaged (upper panels) and peculiar line profiles (lower panels) of the Fe Ka line in the case of MCG -6-30-15 observed by ASCA in 1994 (left 
panels) and 1997 (right panels). Image Credit: Fabian et al., Copyright: PASP, 1 12, 1 145 (2000). 



2.4. Variability of the Fe Ka line 



One of the important features of the Fe Ka line is variabil- 
ity of both, its shape and intensity (see e.g. lBhavani & Na ndra, 
2010, and references therein). Fig. [12] shows one of the most 
illustrative examples for such variations which was found in 
ASCA observations of MCG-6-30-15 during 1994 and 1997 
( Fabi an et all [2000 ) . As it can be seen from Fig. [121 in the 1994 
observation, a very broad profile with a pronounced red wing is 
seen during a period of deep minimum of the light curve (lower 
left panel), compared to the time-averaged line profile shown in 
the upper panel. In contrast, during a sharp flare in the 1997 ob- 
servation, the whole line emission is shifted to energies below 6 
keV, and there is no significant emission at the rest line energy 
of 6.4 keV (lower right panel). Both peculiar line shapes can 
be explained b y large gravitationa l redshift in small radii on the 
accretion disk Fabian et al ] (l2000h . 

Variations of double-peaked emission lines are mainly at- 
tributed to physical changes in the accretion disk, but in case 
of some AGN they cannot be explained by a simple, circular 
accretion disk model, and thus the need for asymmetries in the 
accretion disk arose. Such asymmetrie s are usually introduce d 
in form of a precessing elliptical d isk dEracleous et aL . 1995), 
a circular disk with a spiral arms dLewis et all l2010t) or with 
rotating bright spots (Fl ohic & Eracleousl 12008). 



Iwasawa et al. (2004) detected a variable "red" feature of the 



Fe Ka line at 6.1 keV (in addition to the stable 6.4 keV line 
core) in X-ray spectrum of Seyfert galaxy NGC 3516, observed 
by XMM-Newton satellite (see Fig. [T3~l >. This feature varies 
systematically in the flux at intervals of 25 ks and in energy 
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Figure 13: Light curves of Seyfert galaxy NGC 3516 jlwasawa et ail 120041) 
for: the 0.3 - 10 keV band (top), the Fe Ka line red feature (middle) and the 
6.4 keV line core (bottom). Image Credit: Iwasawa et al., Copyright: MNRAS, 
355, 1073 (2004). 
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Unperturbed line 
Perturbed line 




Figure 14: Left: perturbed emissivity of an accretion disk in Schwarzschild metric IJovanovic & Popovic, 2009). Right: the corresponding perturbed (dashed line) 
and unperturbed (solid line) Fe Key line profiles. Image Credit: Jovanovic & Popovic, Copyright: Black Holes and Galaxy Formation, Nova Science Publishers Inc, 
Hauppauge NY, USA, 249-294 (2009). 



between 5.7 and 6.5 keV. Ilwasawa et al.l (120041) found that the 
spectral evolution of the "red" feature agrees well with hypoth- 
esis of an orbiting spot in the accretion disk. 

Appearances of such bright spots co uld be described by per- 



ing from their accretion disks at relativistic speeds, which im- 
print multiple broad absorption features bl ueward of the emis- 



sion lines (see e.g. IChartas et al. , l2002bh . These additional 



turbations in accretion disk emissivity ( Jovanovic et al. . 2010) 



absorption components, arising from outflowing winds, may 
distort the energy region near the Fe Ka line. Recently, 
Chart as et al. (2009) confirmed the presence of X-ray broad ab- 



and could be caused by several physical mechanisms, such as: 
disk self-gravity, baroclinic vorticity, disk-star collisions, tidal 
disruptions of stars by central black hole and fragme nted spi- 
ral arms of the disk (see e.g. IJovanovic et all l2010l and ref- scale of ~ 10 R g (see also iBrandt et all l2009h . iDone et al 
erences therein). All these phenomena have different occur- 



sorption lines from the quasar APM 08279+5255, varying on 
short timescales of several d ays which implies a s ource size- 



rence frequencies, durations, timescales, strengths, proportions 
and other characteristics. Especially, perturbations of accre- 
tion disk emissivity in form of flares with high amplitudes are 
of great significance because they could provide information 
on accretion physics under extreme conditions. Such flares 
with the highest amplitudes are usually interpreted in terms of 
tidal disruptions of st ars by supermassive black holes (see e.g. 



Komossa et al., 2008, and references therein). Stars approach- 



ing a SMBH will be tidally disrupted once the tidal forces of 
the SMBH exceed the star's self-gravity, and part of the stellar 
debris will be accreted, producing a luminous flare of radia- 
tion which lasts on the timescale of months to years. Although, 
frequency of such events in a typical elliptic al galaxy is very 
low, between 10 -5 and 10 -4 per year (s ee e.g. Jovanovic et all 



2010, and references therein), recently iKomossa et al.l (2008) 
reported the discovery of an X-ray outburst of large amplitude 
in the galaxy SDSS J095209.56+214313.3 which was proba- 
bly caused by the tidal disruption of a star by a supermassive 
black hole. A simulation of such perturbed emissivity of an 
accretion disk in Schwarzschild metric, as well as the corre- 
sponding perturbed and unperturbed Fe Ka line profiles, are 
presented in the left and righ t panels of Fig. [14] respectively 
( Jovanovic & Popovic . 120091) . 

Most AGN contain energetic outows of ionized gas, emerg- 



(2007) found an evidence for a P Cygni profile of the Fe 
Ka line in narrow line Seyfert 1 galaxies and showed that a 
sharp drop in a such profile at ~ 7 keV results from absorp- 
tion/scattering/emission of the iron Ka line in the outflowing 
wind. 

In addition to intrinsic causes such as a disk instability, the 
Fe Ka line variability could also be induced by some exter- 
nal effects, such as ab sorption by X-ray ab s orbers or grav- 
itational microlensing (IJovanovic & Popovic . l2009h . Obser- 
vations of so-called L ow Ionizatio n Broa d Absorption Line 
quasars (e.g Mrk 231 brako et all |2004 and (H 1413+117 



Chartas et al. , 2007|) confirmed the presence of X-ray absorbers 



in these objects. Iwang et aL I (l200ll) detected an absorption line 
at 5.8 keV in nearby Seyfert 1.5 galaxy NGC 4151 and a vari- 
able absorption l ine at the same energy has been discovered 
by iNandra et al. dl999h in NGC 3516. It was interpreted as 
a Fe K resonant absorption line, redshifted either by infalling 
absorbing material or by strong gravity in the vicinity of the 
black hole. A model of the Fe Ka line absorbing/obscuring re- 
gions in form of absorbing mediu m comprised of co ld absorb- 
ing cloudlets was developed by iFuerst & Wul §004). Another 
model, which assumes that absorbing region is composed of 
a number of individual spherical abs orbing clouds, was devel- 
oped by IJovanovic & Popovic d2007h . 
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Figure 15: Left: Simulated images of accretion disk in Schwarzschild metric and the corresponding Fe Ka line profiles. Right: The same, but in Kerr metric for 
maximally rotating black hole, i.e. for a = 0.998. Disk inclination takes the following values from top to bottom panels, respectively: (' = 5°, 30°, 45°, 60°, 85°. 
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Popovic et al.l pOOl) found that a compact gravitational mi- 
crolens could induce noticeable changes in the Fe Ka line 
profile when the Einstein ring radius associated with the mi- 
crolens is comparable to the size of the accretion disk. Also, 
the X-ray continuum could experienc e significant amplifica - 



tion during such micr olensing events (Po povic et al. , 2006b) 



Zakh arov et al. d2004 found that cosmologically distributed 
gravitational microlenses could significantly contribute to the 
X-ray variability of high-redshifted (z > 2) quasars. Indeed, 
microlensing of the Fe Ka line has been reported at least in 
three lensed quasars w ith multiple images: MG J04 14+05 3 4 
dChartas et all l2002ah. OSO 2237+0305 dDaietall | 2003h . 



and H 1413+117 dOshima et all 1200 U IChartas et aU 120041) 



Chart as et al. (2002a) detected an increase of the Fe Ka equiv 



alent width in the image B of MG J0414+0534 which was not 
followed by the continuum and explained this behavior by as- 
sumption that the thermal emission region of the disk and the 
Compton up-scattered emission region of the hard X-ray source 
lie within small er radii than the iron-line reprocessing region. 
Dai et al.l d2003l) also measured amplification of the Fe Ka line 
in component A of QSO 2237+0305. For more detailed dis- 
cussion about observational and theoretical investigations of 
the Fe Ka line var i ability due t o gravitational microlen sing 
see e.g. 



Jovanovicl d2005L 12006): Jovanovic et al.l d2008l) and 



Popovic et alj (2003a b); iPopovic et al.l (l2006al) . as well as ref- 
erences therein. 



the underlying reflected and direct components is crucial in us- 
ing the Fe Ka line as a diagnostic of General Relativity in the 
strong field regime. 

Except the black hole spin, the disk inclination also have sig- 
nificant influence on the line profile, as can be seen from Fig. 
[TBI For small inclinations, i.e. in the case of almost face-on 
disk, the single peak profiles of the Fe Ka line are obtained. 
For inclinations about i ss 30°, the faint "red" peak is the most 
emphasized. It is interesting that thi s value is close to the aver- 
aged value of i = 35°, estimated by Nandra et al.l d 1997b from 
the study of the Fe Ka line profiles of 18 Seyfert 1 galaxies. 

4. Conclusions 

The results of previously mentioned observational and theo- 
retical studies indicate that the broad Fe Ka line which origi- 
nates in vicinity of the supermassive black holes is a powerful 
tool for studying their properties, space-time geometry (metric) 
in their vicinity, their accretion physics, probing the effects of 
their strong gravitational fields, and for testing the certain pre- 
dictions of General Relativity. However, some important sci- 
entific issues still need to be addressed, such as decoupling the 
narrow and broad Fe Ka line components and expl aining their 
different variabilities (see e.g lSulentic et al. , ll998l) . This will 
be one of the crucial steps in determining the real importance 
of relativistic effects on the observed Fe Ka line profiles. 



3. Discussion 

In order to demonstrate the consequences on the Fe Ka line 
of some general relativistic and strong gravitational effects, 
which are significant only near the marginally stable orbit, we 
performed ray-trac ing simulations (using approach proposed by 



Cadez et al. . 1998) of a disk emission from its innermost re- 



gions. For that purpose we analyzed several cases for disk in- 
clination and spin, assuming that the disk extends from the R ms 
to 10 R g , which is the region where the line is most like ly pro- 



duced in most AGN (see e.g Ballantvne & FabianL 2005 ). 

The obtained results are presented in Fig. Q3] As one can see 
from the Fig. [15] for all studied inclinations there are signifi- 
cant differences between the simulated line profiles in the case 
of non-rotating and maximally rotating black hole. In the lat- 
ter case the line profiles are generally wider and more extended 
towards lower observed energies (see also Fig. [TTJ, due to fact 
that the disk extends from the radius of marginally stable or - 
bit, which is much smaller in this case d Jovanovic et all 1201 lb . 
This effect, when observed, gives opportunity to measure black 
hole spin, and thus to have insight into space-time geometry in 
its vicinity. At the same time, as a consequence of gravitational 
redshift, the "blue" peak is more smeared than the "red" one in 
the case of maximally rotating black hole. Therefore, compar- 
isons between the observed and simulated Fe Ka line profiles 
enables us to investigate the effects of strong gravitational field 
predicted by General Relativity. However, since the spectrum 
"reflected" from the disk may contain not only a single iron 
line, but also many fluo rescent lines a nd a continuum reflec- 
tion component (see e.g. iFabian , 2006), accurate modeling of 
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